Strong correlation effects in diatomic molecular electronic devices 
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We present a qualitative model for a fundamental process in molecular electronics: the change in 
conductance upon bond breaking. In our model a diatomic molecule is attached to spin-polarized 
contacts. Employing a Hubbard Hamiltonian, electron interaction is neglected in the contacts and 
explicitly considered in the molecule, enabling us to study the impact of electron interaction on the 
molecular conductance. In the limit where the electron repulsion is strong compared to the binding 
energy (as it becomes the case upon dissociation) electron transmission in strongly suppressed com- 
pared to the non-interacting case. However, the spin-polarized nature of the contacts introduces a 
coupling between the molecular singlet and triplet states. This coupling in turn yields additional 
resonances in the transmission probability that significantly reduce the energetic separation ap- 
pear between resonances. The ramifications of our results on transport experiments performed on 
nanowires with inclusions H2 are discussed. 



PACS numbers: 71.10.-w, 72.10.-d, 72.25.-b 

Molecular electronics [TJ [3 [3 SJ [5] offers new perspec- 
tives for the construction of electronic devices. Even 
though some quantum effects, such as tunneling, are 
sometimes undesirable in conventional electronics, molec- 
ular electronic devices (MEDs) systematically exploit 
quantum effects and are therefore suitable for the con- 
struction of devices on the atomic scale. Another field 
that may contribute to the development of new electronic 
devices is spintronics (for reviews see, e.g., [3E]), where 
spin and charge transport are investigated on a meso- 
scopic scale. Molecular spintronics emerged in experi- 
mental and theoretical studies as a result of recent com- 
bination of molecular electronics and spintronics [SJ [9]. 
Most theoretical approaches to MED's are based on ef- 
fective one-electron theories [10111111131131111151113113 
[T8] and do not properly account for electron correlation 
effects [TTJ [19] . Providing a remedy for this shortcoming 
remains a formidable challenge. Recent work to address 
this issue invokes, for instance, the configuration inter- 
action method to calculate the conductance. While the 
proper inclusion of the infinite contacts remains an un- 
solved problem in this approach, correlation effects are 
accounted for. Non-equilibrium Green's function tech- 
niques have also been employed [20] in conjunction with 
model Hamiltonians to understand correlation effects. 
Here, we follow a similar route. We use the recently 
developed source-sink potential method (SSP) [2TJ[22] to 
a model Hamiltonian involving Hubbard term to address 
the problem of electron correlation in molecular conduc- 
tors. 

One of the most fundamental phenomena undergone by 
molecules is that of bond formation and bond breaking. 
Chemical bonds are sometimes stretched and broken by 
mechanical force in break junction experiments [23} [24]. 
In this letter, we address the breaking of a chemical bond 



in a diatomic molecule, which constitutes part of a MED. 
We provide an answer to the question of how the con- 
ductance of a diatomic molecule changes during bond 
breaking. In a simple molecular orbital picture (appro- 
priate for non-interacting electrons), during the breakage 
of a bond, the binding and antibinding orbital approach 
each other in energy until they become degenerate in the 
separated atom limit. Therefore, electron interaction ef- 
fects become increasingly important on the scale of the 
single-electron energies and a strongly correlated system 
is obtained. This system cannot be properly described 
by the Kohn-Sham reference determinant [25 , 26 , or any 
other single-particle description, appropriate for a non- 
degenerate ground state. As a consequence, the popular 
combination of DFT and non-equilibrium Green's func- 
tion technique [TTJ [27, 28 to model transport in MEDs 
is inadequate. 

The system under consideration is comprised of two 
semi-infinite non-interacting contacts. The boundary 
conditions of the problem are chosen such that in the left 
contact (L) we have an incoming Bloch wave that pro- 
vides an electron current towards the contact-molecule 
interface. At the interface, the incoming wave is par- 
tially reflected from and partially transmitted through 
the molecule. In the present work, we admit only in- 
finitesimal voltage differentials, i.e., we consider the zero- 
bias conductance and assume that the ambient temper- 
ature is zero. We will be concerned with noninteracting 
spin-polarized contacts which contain up-spin electrons 
only. This requirement can be fulfilled experimentally 
by using ferromagnetic leads with parallel configurations 
[9]. Also, external magnets or impurities in the contacts 
can be employed to spin-polarize the contacts [3 HI]- The 
molecule constitutes a subsystem in which the electrons 
repel each other and correlation effects are accounted for 
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by means of a multi-configurational description. By con- 
straining the wave function, we ensure that only a sin- 
gle electron can escape from the molecule into the con- 
tacts. Having two electrons in the contact would mean 
that the wave function in the contact is a multi-electron 
and not a single-electron one. Apart from spin-polarized 
contacts, our model might be relevant to experiments 
[3Q| |3T] , where hydrogen molecules are confined between 
metal wires. The charging energy of the molecule in such 
wires might result in Coulomb blockade and explain the 
appearance of just one open channel. While the con- 
tacts are non-interacting one-electron regions, the inter- 
face between the molecule and the contact requires some 
attention. The interface becomes relevant for configura- 
tions in which one electron is on the molecule and the 
other one is in the contact. Since there are only up-spin 
electrons in the contacts, the electron remaining in the 
molecule has down spin. One might think that the elec- 
trons are either singlet or triplet coupled, but in fact a 
mixture of singlet and triple state is required to account 
for the present electron distribution. This is because the 
contact electron is an up-spin electron which cannot be 
realized with a triplet or a singlet coupled configuration, 
but with an appropriate superposition of the two. This is 
an important point for the basis set selection since in our 
approach, we calculate the stationary states of a model 
Hamiltonian in finite basis set representation. 

The configuratons that we consider as basis functions 
are determined such that they can describe known simple 
limits. For our purposes, the most important such limit 
is the dissociated molecule where the two atoms of the 
molecule are decoupled. The incoming up-spin electron 
in the left contact is completely reflected and it does not 
reach the right contact. Furthermore, due to electron- 
electron repulsion, the up-spin and down-spin electrons 
are localized in the left and right part of the system re- 
spectively. Another obvious limit is the non-interacting 
on. In this case, we have to recover the known transmis- 
sion probability (T(E)) curve (see, e.g. [21 j) exhibiting a 
resonance due to the HOMO and one due to the LUMO 
orbital. One of the key results obtained from our model is 
that the separation between the HOMO and the LUMO 
resonance is drastically reduced as the electron interac- 
tion is introduced. 

We study T(E) of this system as a function of the en- 
ergy with which the electron travels through the contact, 
i.e., as a function of the Fermi energy at zero tempera- 
ture. Alternatively, one might think of a gated device in 
which the Fermi energy of the contacts is kept fixed and 
the molecular part of the MED is shifted in energy by 
a gate voltage. In doing so, we obtain resonances that 
are related to the molecular states. Therefore, before in- 
vestigating the MED, we consider the isolated molecule 
(i.e. sans contacts) subject to electron-electron interac- 
tion. The Hamiltonian of the isolated molecule is given 



by 
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where the operator c^ aa (c aa ) with a=a,b cre- 
ates (destructs) an electron in the molecular sites a and 
b. <j stands for the spin degrees of freedom and n aa is 
the number operator. U is the Hubbard electron-electron 
repulsion parameter and /3m the intramolecular hopping 
parameter. The diagonal elements h a and are set equal 
to zero, thus we will be concerned with a homonuclear 
molecule. We consider a finite Hilbert space in which to 
diagonalize the molecular Hamiltonian. The singlet con- 
figurations that form the basis functions can be divided 
into two ionic ones with both electrons localized on either 
atom a or atom b and a covalent one with one electron 
per atom, i.e., 

*i = -j=\aa>(\U>-\l», 

$2 = |(| ofc > + | 6o >)(|U> - UT>), 

$3 = _L |&&>(|u>-UT>). (2) 

This basis enables us to describe the ground state (dou- 
bly occupied HOMO) and two excited states where one 
electron or two electron are promoted to the LUMO, re- 
spectively. The Hamiltonian in the finite subspace is then 
given by (32] 

/ U V2/3 M \ 
H = V2/3 M V2/3 M . 
V V2/3 M U J 

Solving the eigenvalue equation for the isolated molecule, 
we obtain three eigenvalues, 

Ex = U 

E 2 = ±(u+y/lP + 160* M y (3) 

As U is increased, the ground state raises in energy to 
reach zero for U = oo. In the corresponding wave func- 
tion, there is one electron on each atom, completely elim- 
inating the repulsion energy. Due to this localization of 
the electrons there is no binding interaction between the 
atoms either. Furthermore, since our reservoirs are spin 
polarized, the triplet state of the system 

$ 4 = I(| 06 > - | 6o >)(|U> + UT>)- (4) 

also plays an important role. The energy of this state is 
zero, independent of the value of U. In the limit where 
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Pm goes to zero or where U goes to infinity, the singlet 
ground state and the triplet state become degenerate. 
For the MED this implies that in these limits, a super- 
position of singlet and triplet state determines the trans- 
mission probability through the system. 

We can now turn our attention to the MED, which 
consists of the desribed diatomic molecule attached to 
the left (L) and right (R) spin-polarized semi-infinite con- 
tacts in serial arrangement. The contacts are modeled by 
one-dimensional chains of atoms, nearest neighbor cou- 
pled with strength pc- The diagonal elements of the 
contact atoms are set to zero. L and R in turn, are cou- 
pled to atom a and 6, respectively, through the parameter 
Pcm- Employing the SSP approach [21] to eliminate the 
infinite contacts, we obtain the Hamiltonian 



H = hf +h\ + Vt 



where 
h 



hi = 



and 



E L n LT + T, R n R ^ + Pcmc ] l ^ + H.c. H 
PcMC ] m c h ^ + H.c. + Pmc^c^ + H.c. 

PcMC ] L[ c ai + H.c. + PcMC ] R[ c hi + H.c. 
pMcl i c bi +H.C. 

V U = ^2Un a] n ai . 



(5) 



(6) 
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The operators c^ acr (c acr ) with a=a, b and c^ a 
with /3=L,R create(destruct) an electron in the molecu- 
lar sites a and b and contacts L and R respectively. The 
source El = Pc~ — l+^ e 9 and sink E# = fic e%q poten- 
tial [21] appearing in the Hamiltonian are applied to the 
first atom of L and R respectively. In these potentials, 
q = arccos i^^j- They represent the infinite parts of 

single-electron contacts in a rigorous manner [2T] , 

Building on the discussions above, our finite dimen- 
sional Hilbert space consists of the for configurations 
listed in Eqs.[2]and[4] We additionally include configura- 
tions that are obtained from the described ones by excit- 
ing the up-spin electron from atom a or atom b into the 
contacts. The resulting eight configurations are depicted 
in Fig. [I] The Hamiltonian Eq. [5] is now represented in 
the space of the configurations, 
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FIG. 1: This figure schematically shows the states used to 
construct the Hamiltonian matrix for the full MED. The ar- 
rows represent the spin orientation of the electron and their 
positions depict where the electron is localized within the 
MED. The letters denote whether the diagram corresponds 
to a singlet or triplet. 



Having defined our model Hamiltonian, we can proceed 
to describe how the reflection coefficients are calculated. 
The transmission probability T(e) gives the probability 
of an electron arriving with an energy e to pass through 
the molecule into the right contact. T(e) is related to the 
reflection coefficient r(e) by T(e) = 1 — |r(e)| 2 . r(e) is 
an as of yet unknown parameter in the potential E £ . To 
obtain r(e), we start from the eigenvalue equation 



(9) 



where E is the total energy of the system, e takes a real 
value of our choice. The total energy of the system can 
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,[34], where E 



N-l 



now be written as E = E- 
is the energy of the N — 1-electron system obtained by 
removing the up-spin electron. In the present example 
E = E N ~ X can take on two different values, one being 
the orbital energy of the binding and the other being the 
orbital energy of the anti binding orbital. Since we are 
interested in the state that emerges from the molecular 
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e, where E^ 1 is the 



ground state, we set E = E { 
energy of the binding orbital in the isolated molecule. 
Therefore, Eq. [9] becomes 
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H(e,r)V = {E^~ l + e)*. 
The corresponding secular equation is 

Det(H(e,r)-(E»- 1 +e)) = 



0. 



(10) 



(11) 



Since the value of e represents a boundary condition that 
we specify, Eq. 11 determines the unknown r in terms of 
the variable e. 

First, we verify that our model reproduces the known 
answer of T(e) = 1 the non-interacting, perfect ballis- 
tic conductor. Therefore, we chose a homogeneous linear 
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wire with (3 c = (3 cm = Pm = — 2 and U = 0. Note that 
we use arbitrary energy units as only the relative pro- 
portions of the parameters are of significance. Calcula- 
tion of T(e) according to the above procedure yields the 
desired solution (red dashed curve in Fig. [2]). We also 




FIG. 2: This figure shows T(e) for a homogeneous non- 
interacting wire as a function of Fermi energy e. Two so- 
lutions are found and represented in different colors. The 
second solution (blue curve) corresponds to a wave function 
that exhibits a localized excited state in the " molecule" . 

find a second solution (blue dashed curve) that varies 
strongly. This curve originates from the many-electron 
basis employed for the central region (supermolecule) of 
the MED comprising of the molecule plus the two contact 
atoms explicitly considered. Even when [7=0, we still 
have a many-electron wave function in the supermolecu- 
lar region. This means that we have to expect resonances 
that correspond to local excitations in the supermolecule. 
Such excited states also couple to the contacts and yield 
a T(e) that is neither constant nor equal to one. In the 
localized excited state, the down-spin electron is in the 
anti-bond and the incoming electron passes through the 
HOMO. Since such a local excitation is unlikely to remain 
localized in a system rather than diffusing into the con- 
tacts, we think that the second curve is not observable. 
In our simple treatment, we lack the basis functions in 
the contacts that would permit the local excitation to de- 
cay. While this excitation by itself might not necessarily 
be significant, its impact on the non-excited solution is. 
We pursue this issue for inhomogeneous MEDs, where it 
is realistic to have a strong electron-electron interaction 
in only the molecular part of the MED. 

We now consider an inhomogeneous system, where the 
contact atoms are more tightly bound (3c = —10 than 
the remaining ones ((3cm = —2, (3m = —2). One might 
think of this system as a stretched break junction, with 
a weakened molecule-contact and molecular bond. The 
corresponding T(e) curves for vanishing U are shown in 
Fig. [3j The blue dashed curve on the left, emerging from 
the ground state of the molecule has a resonance at the 
HOMO energy and a second one at the LUMO energy. 
The separation of these orbitals in the isolated molecule 
is 4(2 Pm) and this separation is essentially preserved [22 
in the MED since the value of (3cm is large, this is be- 
cause, we are close to the wide band limit where the 
contacts become featureless. Note that this curve coin- 



cides with the one obtained from a non-interacting the- 
ory (e.g., SSP as described in Ref. [22 j). The red dashed 
curve, emerging from an excited state of the molecule 
shows a resonance at the LUMO energy and a further 
resonance corresponding to a doubly occupied LUMO or- 
bital. The red-dashed curve is absent in an independent 
electron treatment. Even though the interaction param- 
eter is U = 0, we still observe excited states (and the cor- 
responding resonances) for the reasons discussed above. 
As in the example of a homogeneous wire, the excited- 
localized state would probably be quenched by the con- 
tacts and difficult to observe experimentally, unless the 
coupling between the contacts and the molecule is small. 
For finite [7, the excited states couples with the ground 
state and therefore, becomes noticeable in the T(e) curve 
corresponding to the ground state. 




FIG. 3: Left panel shows T(e) for an inhomogeneous wire 
containing a diatomic molecule for [7=0 (dashed curves) and 
for [7=4 (solid curves) as a function of Fermi energy e. Right 
panel depicts the same system for [7=8 

In the interacting system that we are going to consider 
next we set [7 = 4. This choice is motivated by the fact 
that in the isolated hydrogen molecule, [7 has to be twice 
as large as \(3m \ in order to describe the certain properties 
of the hydrogen molecule [35]. For an interacting system, 
with [7 = 4, the ground state is destabilized compared 
to the non-interacting system (see Eq. [3|. Furthermore, 
the separation between the ground- and the first-excited 
state increases to 6.5. Therefore, naively one would ex- 
pect a corresponding increase in the separation between 
the resonances in the T(e) curve. We investigate this 
question with the help of the solid curves in Fig. [3] where 
[7 = 4. Surprisingly, the separation between the first two 
resonances diminishes even compared to the [7 = case. 
This is due to the mixing of states that yield the coincid- 
ing resonances in the T(e) curve of the non-interacting 
case. Furthermore, the second solution (red solid curve) 
is shifted towards higher energies due to the electron- 
electron repulsion. The states in the red solid curve have 
large ionic contributions in which two electrons are local- 
ized on one molecular atom. Increasing [7 even further 
to a value of [7 = 8 yields the curve depicted in the right 
panel of Fig. [3J The gap between the two resonance di- 
minishes further and the resonances in the red curve move 
higher in energy. Eventually, when the two resonances in 
the blue curve start to overlap with each other destructive 
interference sets in, leading to the disappearance of T(e). 
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In Fig. [4] we show the transmission results for U = 32. In 
the large U limit, the ionic states described above acquire 
an infinite energy due to electron repulsion. The two co- 
valent states, having an electron on atom a and a second 
one on atom b in a singlet and triplet spin-coupling, be- 
come degenerate. As a consequence of this degeneracy, 
we get a destructive interference between the molecular 
singlet and triplet state leading to a localization of up- 
spin electron and down-spin electron on atom a or atom 
b respectively. This interpretation can be verified by ex- 
amining the coefficients if the various configurations in 
the wave function. 
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FIG. 4: T(e) for an inhomogeneous wire approaching the 
strong correlation limit (U — 32) as a function of Fermi en- 
ergy e. 

Next, we consider the case where the coupling between 
the molecule and the contacts is weak (3cm < 1- We 
obtain a T(e) curve that resembles the one of Fig. [3j ex- 
cept, of course, the peak widths are reduced. Increasing 
the interaction parameter (U = 4) yields similar results 
to the ones shown in Fig. [3] Again, we observe a dra- 
matic reduction of the separation between the first two 
resonances as compared to the non-interacting case. 

Now we turn to the dissociation limit, where the cou- 
pling between the molecular atoms (/3m) approaches zero. 
In this case, on the scale of the (3m, the importance of 
the repulsion parameter U increases and we enter the 
strong-correlation regime. Therefore, we expect the same 
qualitative behavior as in the strongly interacting case. 
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FIG. 5: The left plot shows T(e) for weakly coupled ((3 m = 
— 1) molecular atoms. The dashed curves correspond to non- 
interacting system. The solid curves are obtained for the same 
system with U — 4. The right plots illustrates T(e) for very 
weakly coupled molecular atoms ((3 m — —1/2) with [7 = 4. 

In the case where U = and (3m = — 1 (Fig. |5| all 
the resonances are compressed into a smaller energy in- 
terval as expected. This implies that electron interaction 



is going to have a large effect, since all these close res- 
onances are going to be coupled. The solid blue curve, 
where U = 4 yields two coinciding resonances, or equiv- 
alently, the gap between the first two resonances dis- 
appeared completely. For cases where the Fermi level 
is located between the HOMO and LUMO level of the 
non-interacting case, a strong conductance enhancement 
would be obtained. As already explained, this enhance- 
ment is due to the indirect impact of the molecular triplet 
state. Reducing the absolute value of /3m even further 
((3m = — 1/2) in the right plot of Fig. [5]) leads to a strong 
suppression of the T in the interacting case as opposed to 
the non-interacting case. We are again approaching the 
strongly-correlated limit and naturally, the T(e) curve 
resembles the one of Fig. [4j The preceding discussion 
underlines the statement made in the introduction, that 
bond breaking turns the system into a strongly correlated 
one. Therefore, even if in reality, the value of the parame- 
ter U is fixed, stretching the bond amounts to an effective 
increase of the parameter relative to the atom-atom cou- 
pling. In general, one would expect that electron-electron 
interaction effects reduce T(e) [19 . This is indeed what 
we observe, except in cases where the electron repulsion 
mixes in further states as it is the case here with the 
molecular triplet states. 

In conclusion, our understanding of correlation effects 
on the transmission probability of molecular conductors 
is still quite limited. In this letter, we provided a qual- 
itative model of correlation effects. We focused on the 
molecular part of the MED and considered the simplest 
possible way of introducing electron interaction effects 
therein. In the interacting case, we uncovered a new, un- 
expected resonance that originates from the lowest triplet 
state of the isolated molecule. In the presence of electron 
repulsion, this state is coupled to the molecular ground 
state to drastically reduce the gap between the HOMO 
and the LUMO resonance that would be observed in a 
non-interacting system. 

In our model, the contacts were assumed to be spin 
polarized. Nonetheless, we argue that our model calcula- 
tions pertain to realistic systems. Indeed, recent experi- 
ments [30j [31] in which H2 molecules are inserted between 
gold wires give rise to the speculation that there is only 
one spin channel open. Spin-polarization in the contacts 
appears to be attainable with experimental techniques as 
well. The fact that the conductance is suppressed grad- 
ually in our calculations in the strongly correlated limit 
as a function of (3m instead of collapsing abruptly serves 
as further evidence that our results are in line with the 
experimental data [30l |36l |37l [38] - Therefore, we have 
reason to believe that our model sheds light on electron 
interaction effects despite its seemingly unsophisticated 
nature. 

We gratefully acknowledge financial support provided 
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